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Spatial bandlimitedness of scattered 
electromagnetic fields 

Uday K Khankhqje and Kushal Shah 


Abstract —In this tutorial paper, we consider the problem 
of electromagnetic scattering hy a hounded two-dimensional 
dielectric object, and discuss certain interesting properties of 
the scattered held. Using the electric held integral equation, 
along with the techniques of Fourier theory and the properties 
of Bessel functions, we show analytically and numerically, that 
in the case of transverse electric polarization, the scattered fields 
are spatially bandlimlted. Further, we derive an upper bound 
on the number of incidence angles that are useful as constraints 
in an inverse problem setting (determining permittivity given 
measurements of the scattered field). We also show that the 
above results are Independent of the dielectric properties of the 
scattering object and depend only on geometry. Though these 
results have previously been derived in the literature using the 
framework of functional analysis, our approach is conceptually 
far easier. Implications of these results on the inverse problem 
are also discussed. 

Index Terms —Electromagnetic scattering by nonhomogeneous 
media. Inverse problems. Integral equations. 


I. Introduction 

Scattering is a phenomenon where some form of a travelling 
wave excitation (light, sound, etc.) deviates from its original 
trajectory due to a change in the properties of the medium 
along its path IT]. In the context of electromagnetic waves, 
given a field incident on an object of known permittivity, it is 
quite straightforward to calculate the scattered field in various 
directions. In cases where this cannot be done analytically, sev¬ 
eral computational methods can be employed. This is known as 
the forward scattering problem. The inverse problem, however, 
is more challenging. It consists of determining the unknown 
spatial permittivity of an object based on measurements of the 
scattered field Q. 

In order to understand the properties of the scattered field, 
Bucci et al. El, El considered the electric field integral 
equation (EFIE) approach, and noted that the integral operator 
in this case is compact. By invoking a theorem due to 
Kolmogorov and Fomine 0 concerning the properties of such 
an operator, it was deduced that the scattered field has a finite 
dimensional representation. Further, the singular values of the 
operator rapidly decay after a certain threshold, a property 
attributed to the analyticity of the operator 0. Thus, it was 
concluded El that the scattered field can be represented by 
a finite number of singular vectors, each associated with a 
singular value. This critical number was referred to as the 
degrees of freedom of the scattered field, and in the two- 
dimensional case of a circular observation domain bounding 
a circular scatterer of radius a, this number was found to be 
equal to 2ka, where k is the free space wavevector 0. Also 
see I 2 I for a lucid derivation of this decomposition. 

U. K. Khankhoje and K. Shah are with the Department of Electrical 
Engineering, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 
110016, India. E-mails: uday@alumni.caltech.edu and kkshah@ee.iitd.ac.in, 
respectively. 


The mathematical machinery used in the works of Bucci 
et. al is formidable, as it is rigorous. Instead, we present a 
much simpler route to the same results by invoking certain 
elegant properties of Bessel functions. We start in Section 
[II| by discretizing the EFIE, and derive expressions for the 
discrete Fourier transform of the field scattered by a bounded 
dielectric object in the case of transverse electric polarization. 
Using a key property of Bessel functions of fixed argument 
(demonstrated in Appendix |^, namely that their amplitude 
monotonically decays to zero as the order is increased beyond 
a threshold value, we present our main results on the bandlim- 
ited nature of the scattered fields in Section [111] In this section, 
we also present numerical validation with the finite element 
method, and consider the case of an object illuminated by 
multiple incidence angles. We conclude in Section IV with a 
discussion on the hardware implications of our results, and 
make a few other connections. 


Observation m 



Fig. 1. A dielectric object T) bounded within a circle of radius a is illuminated 
by an incident field. The scattered field is measured on the circle E of radius 
po- refer to discretized source and observation points, respectively. 


II. Methods 

A. Deriving the scattered electromagnetic fields 

Consider a dielectric object bounded in a two-dimensional 
region of space, D, whose relative permittivity, er{f) is a 
function of space, immersed in a medium of constant (and 
real) relative permittivity, e^. In this situation, if we consider 
a transverse electric (TE) polarization {H^ = 0), the z- 
component of the electric field obeys the Helmholtz equation; 

^^E,{f) + eer{f)E,{f)=Q (1) 

where k = 27r/A is the magnitude of the wavevector in free 
space, and er{f) in the relative permittivity at r. The above 
equation can be recast as an integral equation in terms of the 
Green’s function for the homogeneous medium (with uniform 
permittivity, Cb), G(r, F') = [—j/A)E[lf^{kb{f—r')), where 
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( 2 ) 

Hq is the Hankel function of second kind and zero order, 
and kf) = k^/cb, as; 


= Elif) + kl ^)x{nEA^) dT^ ( 2 ) 

V 


where E], is the incident held, and x(f^ = ej.(f)/^ is 
referred to as the dielectric contrast. The above EFIE is a 
Ereedholm integral equation of the second kind. 

Assume that the dielectric object, T), is surrounded by 
a concentric observation circle, E, of radius po shown in 
Figure To compute the scattered held {E^ — El) on T, we 
discretize Eq. 0 by dividing the region V into N equi-sized 
cells of uniform dielectric constant and T into M equispaced 
points. Following well known techniques for solving such 
equations n, a, the scattered electric held (v) at the to* 
observation point on F is related to the total held (u) at the 
N points of V: 


JV-l 

^ ^ {kbPnin)Xn'^n^ where (3) 

n—0 

TT 

go =-jhan^Jiihan) (4) 

where is the distance between the to* (observation) point 
on r and n* (source) point in T>, an is the radius of the circle 
with the same area as the n* cell of 22 IH with dielectric 
contrast Xn, and Ji is the Bessel function of the hrst kind and 
hrst order. Note that without loss of generality, oq = ai = 
• • • = a^v-i and po does not depend on n. 

We now approximate the scattered held in a far-held setting 
(po > A). The cosine rule gives us that p^n = {^n + 
Pq — 2r„po cos 9mnyG, where 9mn is the angle between the 
position vectors corresponding to the source and observation 
points, fvi and fm, respectively, and |fv„| = po- Thus, in the 
far held, at least to hrst order, pmn can be approximated as 
Pmn ~ Po in the amplitude term and as Pmn ~ Po—fn cos 9mn 
in the phase term, following which (also using the large 
argument approximation of the Hankel function, Hq \z) ~ 
\/2/(7rz) exp (- j {z - 7r/4))): 


H'^\kbPmn) 


■ exp 


irkbPo 
exp {jkbrn cos 9m, 
ho exp {jkbrn cos 9mn) 


? - hpo 


(5) 


where ho is a constant independent of to when the observation 
points are on a circle (the terms in the curly bracket above). 
Finally, the far held approximation for the scattered held can 
be written as; 

N-l 

goho XnUn exp {jkbTn COS 9 mn ) (6) 

n—Q 


B. Discrete Fourier Transform of the Scattered Fields 

We are now in a position to consider the M-point discrete 
Fourier transform (DFT), v G C^, of the scattered electric 


held, V, as obtained in Eq. The /c* Fourier component is 
given by; 

M-l 

m—O 

N-l 

goho exp (jkbrn cos ^mn )Xn'Un' (7) 

n—0 

A simple reordering of the order of summations in 0 
reveals that the inner summation is, in effect, the DFT of a 
plane wave sampled on a circle of radius r^; 

N-l 

Vk = goho'^XnUnX- 

n—0 

M-l 

exp{-j2TTk^^ exp{jkbrnCOs9mn)- ( 8 ) 

m—0 

This is so because the inner summation, in which the source 
position, fn, is hxed, contains 9mn in the exponent. This 
9mn can be expanded as 9mn = {2TTm/M — 9n), where 
9n is the angular position of the n* source point. Since the 
measurement position, r^, goes around a circle of radius po 
as TO goes from 0 to (M — 1), it is evident that 9mn evenly 
samples points spanning 27r. 

III. Results 

A. Bandwidth: Dependence on size 

The DFT of the scattered held as derived in Eq. ([^ can be 
simplihed using the Jacobi-Anger expansion, as follows (see 
Appendix for details) ; 

N— 1 CXD 

Vk = gohoM y] XnUn j^-1^ X 

n—0 q— — oo 

exp [-j{k - qM)9n]Jk-qM{kbrn)- (9) 

A careful examination of the inner summand of Eq. 0 shows 
that for a hxed point in region D (i.e. hxed n), the argument 
of the Bessel function is constant (equal to kbCn) and only its 
order (equal to k — qM) changes with q. 

The Bessel function, Jk{x), has the property that Jk{x) ~ 0 
for |fc| > |"2|a;|] (see Appendix 0- Further, since the maxi¬ 
mum possible value of r„ is a (for a dielectric bounded within 
a cylindrical region of radius a), it naturally follows that for 
k > \2kba~\ the Fourier component, Vk will be negligible. 
In other words, the DFT, v, is bandlimited to this value 
r] = \2kba~\. It must be mentioned that this bandlimit cor¬ 
responds to an effective bandwidth, as the Fourier coefficients 
for p < fc < M — p are negligible, but not identically zero 

n. 

We thus arrive at the same result as Bucci et. al 0 regarding 
the degrees of freedom of the scattered held in terms of an 
effective bandwidth. We note that our approach of discretizing 
the EFIE is similar to one previously proposed ifTOll . wherein, 
starting from a series representation of the Green’s function, 
a truncated Fourier series for the scattered electric held is 
obtained; we essentially extend this idea and estimate the 
truncation number in terms of p. 
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B. Bandwidth: Dependence on permittivity 

Since the internal field coefficients, and the contrast, 
Xn, are upper bounded in magnitude, the Fourier components 
in Eq. that are zero, will continue to be zero, regardless 
of the particular values of the field and contrast. 

To see this quantitatively, consider the fc* Fourier coefficient 
of the scattered field from Eq. We consider the first half 
of the total M coefficients; this is sufficient since the DFT is 
symmetric (see Appendix |B-A[ ): 

N-l 

Vk = gohoMj’^'^ iXn)[unexp{-jk9n)Jk{hrn)] (10) 

n—0 

Assuming the contrast to be bounded such that \xn\ < Vn, 
and the field to be bounded such that |u„| <7, Vn, 
we apply the Cauchy-Schwartz inequality (IX^On^nP < 

n 

^|a„p^|5„p) to the above relation to obtain 

n n 


N-l N-l 

\vk\'^ < \gohoM\'^'^ \xnf \unexp{-jk9n)Jk{kbrn)\^ 

n—0 n—0 

N-l 

< \gohoM\'^N{KN'^ \Jk{kbrn)\'^ (11) 

n—0 

Thus, if the magnitude of the order of the Bessel function, 
|fc|, is large enough such that Jk{kba) ~ 0 , we can see 
that Ififcl ~ 0 as long as are finite. In other words, the 
bandlimited nature of scattered field is independent of the 
object permittivity, and this bandlimit depends only on the 
size of the object relative to wavelength, kbO. 


C. Numerical Results 


We use a two dimensional vector-element based finite 
element method (FEM) im to compute the scattered electro¬ 
magnetic fields in two different configurations. In both con¬ 
figurations, the scattering object is confined within a cylinder 
of radius r = 2A, the fields are computed on a radius of 
r = 3A, the computational domain is terminated by applying 
a radiation boundary condition at a radius r = 4A, and the 
TE-polarized incident field makes an angle of 0° with the +x 
axis. For the numerical convergence of the FEM solution, the 
domain discretization must be on the order of A/20, which 
results in a mesh having approximately 90,000 elements. In 
the first configuration, the cylinder has a uniform permittivity, 
Cr = 4, while in the second configuration, we allow the 
permittivity of the cylinder to be random such that for each 
element is a uniform random variable from 1 to 10 . 

The scattered fields in both cases are shown in Figure 
and their corresponding spatial DFTs in Figure We 
find that regardless of the constitutive permittivity of the 
scattering object, the DFT is bandlimited, and that the band 
limit matches very well with the analytical prediction of 
77 = [2A:t,a] from Section [III- A Here, a = 2A , which gives 
77 = 26 (see Figure |^. It is interesting to note that while 
the observation circle (r = 3A) is not in the far-field of the 
scattering object (a = 2A), the predicted cut-off matches very 


well with the FEM results (which do not assume any far- 
field approximation). This is because the idea of an effective 
bandwidth as derived by Bucci et al. El does not require a far- 
field approximation, even though we assume it here to simplify 
the analysis. 



Fig. 2. FEM-computed real pai't of the scattered electric fields as a function 
of scattering angle as observed on a circle of radius r = 3A for a constant 
permittivity (€r = 4) cylinder (blue curve) and for a cylinder with random 
permittivity (red curve). In the latter, the permittivity of each element in the 
cylinder is a uniform random variable such that er G [1,10] (see top-left 
inset of Figure]^. 
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Fig. 3. Absolute value of the DFT of the scattered fields as shown in Figure 
[^for the two cylinders with constant (blue) and random (red) permittivity. 
The top-right inset shows all Fourier coefficients, clearly demonstrating the 
symmetric and bandlimited nature of the fields for both cylinders. The main 
graph shows the first 50 Fourier coefficients, and it is evident that beyond 
k = 26 the Fourier coefficients have negligible magnitude. The inset on the 
top-left shows the permittivity of the random cylinder as a function of space. 


D. The case of multiple incidence angles 

So far, we have considered the dielectric object to be 
illuminated by a single electromagnetic wave, and found that 
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the scattered field can be represented in terms of 2kba Fourier 
components. In other words, in an inverse problem type 
setting, there will be 2ki,a linearly independent relations that 
can be used to relate the known value of the scattered field to 
the unknown permittivity of the object. In the language of the 
formalism presented in this paper, this amounts to determining 
Xn, n G [0, 1] given the measurements Vm, rn G [0, M—1] 

from Eq. Q. 

A natural question that arises is: by what method can more 
information be gained about the object? It is conceivable that 
by introducing multiple incidence angles, further constraints 
can be imposed on the unknown object permittivity. However, 
we show that the number of such constraints can not be 
increased linearly with the number of incidence angles, and 
that there is an upper limit on the number of useful angles. 

To see this, we go back to the integral equation for¬ 
mulation as in Eq. and note that this equation, when 
discretized, leads to a well-determined system of equations 
0, denoted by F u'’ = F, where F is a “forward” matrix, 
consists of the internal field values for the incidence 
angle, and the incidence angle appears in the known vector, 
F. Explicitly, at the fc* source point in the object domain, 
6^ = exp[—jkbPk cos6ik], where pk is the distance of the 
source from the origin and 9ik is the angle between the source 
position vector, fk and the incident wavevector, k^. We assume 
that the angular spacing between multiple incidence angles is 
uniform. The above matrix system has a solution, giving the 
field at the n* source point as; 


N-l 

= E[^ 


-11 


J nk 


exp[-jkbPk cos e^k] 


( 12 ) 


Noting the similarities between Eq. (j^ and Eq. (12i, it is 
at once clear that if the DET of the above expression is taken 
w.r.t. the incident field index, i, the result would be a band 
limited expression, just as was shown with the DET of Eq. 
w.r.t. the observation point index, m in Section II-B Extending 
the analogy further, it is seen that this band limit is given by 
the expression; 


rjp = m■^ix\2kbPk^ = \2.kbd\. (13) 

Pk 

Thus, only upto rjp incidence angles are useful in imposing 
independent constraints on the unknown permittivity of the 
object. Beyond this number, the field can be reconstructed 
using rjp sampled values and a suitable interpolation scheme 
lEl; no new information can be gained. We note that this 
number depends only on the relative object size, and does not 
depend on the object permittivity. We also note that this result 
has previously been derived in the framework of functional 
analysis by Bucci and Isernia 0. 


IV. Discussion 

A. Implications on the Forward scattering problem 

The Nyquist-Shannon sampling theorem m states that 
for a function that is bandlimited to a maximum frequency 
component, fb, the minimum sampling frequency required to 
reconstruct the function is given by fs = 2fb. Applying this 


theorem to the scattered electric field, which is known to be 
bandlimited, implies that for a fixed object size, a, it is optimal 
to make AkbU equally spaced measurements of the scattered 
field. In an experimental setup it may be necessary to measure 
the fields scattered due to an object. It is common to have a 
dedicated antenna for each measurement; we thus have a lower 
bound on the hardware complexity of the experiment. Eurther, 
if it is desired to estimate the scattered field at more points 
than these, a suitable interpolation scheme can be applied to 
the sampled values of the field ifl^ . 

A setup such as that described above—surrounding an 
object with several antennae—is fairly common in experi¬ 
ments which perform breast cancer detection using microwave 
imaging techniques El. Although our results apply to a 
two-dimensional geometry, the same approach can be easily 
applied to the three-dimensional problem. 


B. Implications on the Inverse problem 

In the inverse scattering problem, the measurements of 
the scattered fields are typically noise corrupted. Also, it 
has been shown that the scattered field is bandlimited to 
2kba. Thus, a robust strategy for determining the unknown 
permittivity of an object would be to: (i) Take the DET of 
the measurements (i.e. obtain Vk, k G [0,M — 1]), (ii) apply 
a bandpass filter which zeroes out all Vk components for 
\2kba~\ < k < M — \2kba~\, and (iii) use the remaining 
2\2kba \ components of Vk in Eq. (j^ for the purpose of solving 
the inverse problem, i.e. determining the dielectric contrast, 
Xn- More constraints on Xn via Eq. Q can be imposed by 
increasing the number of incidence angles, but as shown earlier 
(Section |III-D[ ), the maximum number of useful incidence 
angles is limited to \2kba] as well. 

Thus, the property of bandlimitedness can be used to make 
the inverse problem noise-resilient to a certain degree. 

Appendix A 

Bessel Eunctions of Eixed Argument 

In this section, we consider the properties of the n* order 
Bessel function for a fixed argument, x, i.e. Jn{x), where 
n, X G K. We are primarily interested in explaining a key 
feature of Eigure the observation that J„(x) has a steep 
and monotonic decay in amplitude beyond a certain threshold 
value of n for a fixed x. This threshold value is subsequently 
identified. 

We start by considering the integral representation of J„ (x) 
for integer n as 

1 r 

Jn{x) = — / cos(nr — xsinr) dr. (14) 

Jo 

The argument of the cosine term in the above expression, 
^(r) = nT—x sin r has been plotted in Eigure|^for two values 
of n at X = 27r; n = 6 and n = 12. It is clear from Eigure 
that while J6(27r) is far from being zero, Ji2(27r) « 0. This 
is graphically illustrated in Eigurewhere cos())(t) is plotted 
for both values of n. It is evident that the area under the curve 
is clearly non-zero in the case of n = 6 (bottom right of the 
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Fig. 4. Bessel function Jn{x) as a function of n for x = 27r. 



Fig. 5. Phase 0(r) = nr — atsinr as a function of t for a; = 27r and two 
values of n: n = 12 and n = 6 (left), cos(/>(r) for n = 12 (top right), and 
cos </>(r) for n = 6 (bottom right). 


figure), where as the positive and negative areas seem much 
better balanced in the case of n = 12 (top right of the figure). 

The above argument can be presented in a more abstract 
manner, as follows. For a fixed x, and as n increases, the 
value of \n/x\ exceeds 1 and grows. The phase function, 
^(t) = n{T — ^sinr) has a linear component and a sinu¬ 
soidal component, and as \n/x\ grows, the linear component 
dominates. As a result, the intervals in r (denoted by Ari) 
where 4){t) goes from tott to tott + begin to approach the 
length of the intervals in r (denoted by At 2 ) where (j){T) goes 
from m7r+ ^ to (m+ l)7r. The area under the cos (j){T) curve 
for the above two intervals, thus become approximately equal 
while being opposite in sign, i.e. for m G [0, n — 1] : 

/ cos (/)(t) dr « — / cos 4 >{t) dr. (15) 

J Ari J At2 

Thus by breaking up the total integral of Eq. into such 
intervals, we arrive at the result that as \n/x\ grows beyond 
1, Jn{x) monotonically approaches zero. By noting that this 
behaviour depends simply upon |n/a;|, and by inspection of 
Figure 1^ it can be said that Jn{x) « 0 for \n\ > [2|a;|]. 
The modulus on n come from the observation that | J„(a;)| = 


\J-n{x)\. 


Appendix B 

Discrete Fourier transform of a plane wave 


Consider a z-polarized plane wave propagating in the x — y 
plane, measured on a circular contour in the plane. Let there 
be M observation points on this circle of radius rg, giving an 
observed electric field vector, x, as 


x[m\ = exp 


-jfcgro cos 



m e [0,M- 1], (16) 


for a plane wave travelling along the +x axis. The dis¬ 
crete Fourier transform (DFT) of this vector is given by 

M-l 

x[k]= ^ exp[—)/27rfcm/M]x[TO], for fc S [0, M—1], whose 

m—0 

properties are now considered. 


A. Symmetry of the DFT 

We can simplify it’s symmetric counterpart, x[M — k] by 
noting that for the given x, we have x\m\ = x\M — m], 

M-l 

which leads to: x[M — k] = ^ exp[j27rA:m/M]a:[M — m]. 

m—0 

Substituting M — m = ^, we get 


x[M — k] = ^ exp 
1 

H exp 


fi=M 


—j2TTklI 

M 


j2'Kk{M — /i) 

M 

= x[k] 


Af] X 


(17) 


The last equality follows because the summand at /i = M 
is the the same as the value for p = 0 since the field x is 
periodic with period M. 


B. Band Limited Nature 

To arrive at the band limited nature of the incident field, we 
use the Jacobi-Anger expansion, which states that 


OO 

exp(jzcos6») = ^ fjp{z)exp{jp9), (18) 

P— — 00 


and the following expression for the discrete delta function. 


d[k - k'] 


1 

M 


M-l 

m—O 


j2'Kmk 

M 


exp 


j2'Kmk' 

M 


(19) 


where {k,k'} G [0, M — 1]. Observe that the transformation 
k' s- k' + qM, where q G h, leaves the above equation 
unchanged. 
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The fc* Fourier component of x from Eq. dl 
applying Eqs. (18i and ( [T^ , is 

M-l 

x[k] = exp 


after 


n=0 


—j2'Kmk 

M 


OO 

X X! J^'4)(“*o?’o)exp 

P— — 00 
OO 

= X! fJp{-^o'^o)M5[k-{p + qM)\ 

P— — 00 


j2'Kmp 

M 


OO 

= ^ f-^^MJk-qM{-koro) ( 20 ) 

q——oo 


In other words, the fc* Eourier component of x comprises the 
/c* order Bessel function and it’s M—shifted orders, all of the 
same (fixed) argument. 

Erom the analysis of the properties of Jfc(—fegro) as a 
function of k in Appendix it is clear that the DET of 
the incident plane wave is band-limited. This is because for 
k > [2A:oro], Jfc(—fco^o) ~ 0, and therefore x\k] ^ 0. Thus 
the incident field on a contour of radius rp can be represented 
by \2korQ] coefficients in the DET basis. 
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